LU factorization method in matlab